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Abstract. We present the results of NTT/VLT UBV imaging of a 260 arcmin 2 region containing the Galactic Luminous Blue 
Variable WRA 751, in search for its birth-cluster, i.e. a cluster of young and massive stars spatially and physically associated 
with it. On the basis of the classical reddening-free parameter Q, we have identified a sample of 24 early-type stars with 
colours typical of spectral types earlier than B3. Interestingly, these stars are clustered within a radius of 1' from WRA 751, 
corresponding to about 1% of the imaged field. These stars tightly distribute around (B — V) - 1.67, which in turn defines a 
mean extinction A v ^ 6.1 mag. The 5 brighter (V < 16.2) and bluer (Q < -0.9) stars of the sample have been subsequently 
observed with FORS1 and classified as 3 late O- and 2 early B- stars. The absence of stars earlier than 08 indicates an age 
of the cluster older than 4 Myr, although it could be due to an incomplete sampling of the upper end of the main sequence. 
Nevertheless, the detection of OB stars of class I certainly indicates an age of a few million years. At an assumed distance of 6 
kpc, we estimate a cluster radius of 3.4 pc and a total mass of 2.2 x 10 3 M G . Our discovery is only the second known instance 
of a Galactic Luminous Blue Variable associated with its birth-cluster. 

Key words. Stars: early-type, evolution, fundamental parameters, Hertzsprung-Russell (HR) and C-M diagrams 



petition among stellar gravity, radiation pressure and rotation 
Massive stars (M, > 40 - 60 M , i.e. progenitors of super- (Maeder & Meynet 2 Q00). Unfortunately, observations have 
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k*" 1. Introduction Theoretical predictions indicate that mass loss is a com- 

% . . 

Ctf ■ novae of type II) play a key role in galaxy formation and failed so far to quantitative i y constrain the rotational velocity 
evolution. Their stellar winds and end-of-life explosions set of evolved massive stars Mass loss becomes extreme in me 
the kinematics and chemistry of the interstellar medium, and post . main sequen ce evolution, when stars reduce their mass by 
power galactic fontains and outflows, which, in turn, signif- a factor of > 4 and approach the final stage of SU p ernoV a. It is 
icantly affect galaxy evolution (i.e. feedback, Dekel & Silk, believed that massive stars remove me ma j ority of th eir enve- 
1986). This is true at any lookback time; however, beyond the lope during the Luminous Blue Variable (LBV) phase, via giant 
Local Group, individual stars are no longer resolved and one eruptions that produce circumstellar, chemically-enriched neb- 
only observes their effect integrated over the entire galaxy. To ulae (i e Eta Carinae , Humphreys & Davidson 1994). 
interpret these observations, population synthesis codes are re- 
quired (i.e Bruzual & Chariot 2003, Leitherer et al. 1999). The 

accuracy of these models relies on our detailed knowledge of Nebular expansion velocities and nebular sizes set the life- 

the first step in the ladder, the massive stars evolution, and, in time of the LBV phase to ^ 10 4 yrs (Nota 1997), and this is the 

particular, the dependence of mass loss on initial stellar mass reason why only about 30 candidate and confirmed LBV stars 

and time (cf. Langer et al. 1994). This is still under debate. are known in the entire Local Group (Humphreys & Davidson 

1994, Nota 1997). Of these, eleven belong to the Milky Way 

Send offprint requests to: A. Pasquali and onl Y one confirmed LBV is located in a star cluster (Eta 

* Based on observations obtained with ESO telescopes on La Silla Carinae in Trumpler 14/16, Massey et al. 2001), and a candi- 

and Paranal, Chile, under programmes 271.D-5008, 71.C-0432, and date LBV candidate, VI Cygn No. 12 is associated with Cgynus 

073.D-0037 OB2. In the LMC, only one confirmed LBV out of 10 is a mem- 
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ber of a cluster of massive stars (S Doradus in LH 41, Massey 
et al. 2000). 

The presence of a cluster of OB stars associated with an 
LBV star gives us the possibility to infer the distance, age 
and mass of the LBV progenitor by fitting the cluster Zero- 
Age Main sequence and its turn-off. In this way, Massey et al. 
(2001) have derived a progenitor mass of about 120 M and an 
age of about 1 Myr for both LBV stars. The authors have thus 
suggested that LBVs are a normal stage in the evolution of the 
most massive stars. 

To further test this hypothesis better statistics are required, 
and the discovery of the birth-cluster of a LBV star is thus 
extremely valuable. In this paper we report on the discovery 
of a second case of association between a confirmed galac- 
tic LBV and an OB cluster. The LBV star is WRA 751 in 
the Galactic Carina arm. Long-slit spectra of its circumstel- 
lar nebula (Nota, private communication) revealed the pres- 
ence of three emission-line (i.e. Ha) stars in the close vicin- 
ity of WRA 75 1 . Unfortunately, the spectra only covered the 
wavelength range 6200 - 8000 A, and the 4000 - 5000 A inter- 
val, needed to derive reliable spectral types, was missing. We 
have therefore followed-up these early observations with UBV 
imaging of a region centred on WRA 75 1 to identify other pos- 
sible young, massive stars on the basis of their U BV photome- 
try. We do confirm the existence of moderately reddened early- 
type stars and present spectra of a selection of them, which we 
classify as late O- and early B-type stars of various luminosity 
classes. The absence of stars earlier than 08 indicates a cluster 
age older than 4 Myr. 

The observations and data reduction are presented in Sect. 2, 
while the identification of the WRA 751 cluster is described 
in Sect. 3. The photometric and spectroscopic properties of the 
cluster stars are discussed in Sect.4, and conclusions follow in 
Sect.5. 

2. Observations 

2.1. Imaging 

Our first set of imaging observations was obtained on the night 
of 1/2 April 2003 using EMMI, the visible imager and spec- 
trograph at the ESO New Technology Telescope on La Silla. 
The observations were designed so as to include the region 
immediate to WRA 751 as well as an extended, diffuse neb- 
ulosity visible in the UK Schmidt Telescope Ha survey (Parker 
1997) extending to its Southwest. An area of 18'. 6 x IV A cen- 
tered on the coordinates a(2000) = ll ft 08 m 15 s 5, 5(2000) = 
-60°47'35" was covered in a mosaic of pointings through the 
U, B, and V filters. The U and B frames were obtained on 
the EMMI blue arm, which uses a 1024 x 1024 pix 2 detector 
yielding a field of view of 6'. 2 x 6'. 2 per frame with a scale of 
0".37 pix -1 . The V frames were obtained with the EMMI red 
arm, which uses an array of two 2048 x 4096 pix 2 detectors 
yielding a field of view of 9'.1 x 9'. 9. We used 2x2 binning 
in the red arm, resulting in a scale of 0".333 pix 1 . The field 
was covered in a sequence of ten pointings in each filter con- 
sisting of telescope offsets of 3' westwards, followed by off- 
sets of 5' alternatively northwards and southwards. The offsets 



ensured a broad overlap between adjacent frames resulting in 
each sky position being recorded by at least two (U and B fil- 
ters) or three (V filter) pixels. Exposure times for each frame 
were 120 s (t/), 45 s (B), and 10 s (V). Additional observa- 
tions were obtained of the Landolt (1992) standard field around 
PG 1047 +003 in order to determine the photometric zeropoints 
in each filter. Characteristic extinction coefficients for La Silla 
were taken from the observatory Web pages, and photomet- 
ric system transformation coefficients were obtained from the 
EMMI User Manual. 

The NTT observations were reduced using standard IRAF 1 
tasks to carry out the usual steps of bias subtraction, flat-field 
correction using twilight sky flat field frames, removal of cos- 
mic ray hits, and construction of a bad pixel mask based on 
the pixel value histograms of bias and flat field frames. Precise 
offsets among the images composing the mosaics were deter- 
mined using the position of the numerous stars in the overlap 
areas between consecutive frames as a reference. The images 
were then combined into a single mosaic of the field by av- 
eraging the pixel values at each sky position, discarding those 
flagged in the bad pixel mask. Finally, the mosaics obtained in 
each filter were registered and cropped to the common area. 

A second set of observations was obtained using FORS2, 
one of the visible imagers and low-resolution spectrograph in- 
struments at the ESO Very Large Telescope, in Service Mode 
on the night of 7/8 March 2004. The instrument detector is 
composed by mosaic of two 4096 x 2048 pix 2 detectors yield- 
ing a pixel scale of 0".253 pix -1 with the 2x2 binning 
used. A dither pattern focusing on the northern region of the 
WRA 751 cluster, which was not included in our previous 
NTT observations, was observed through the UBV filters. We 
used the pipeline-reduced individual images included by the 
Observatory in the data package (Silva & Peron 2004) and 
stacked them into a single frame per filter for a given expo- 
sure time. Two exposure times per filter were used: 5s and 
120s (U), and Is and 60s (B, V), so as to increase the dynamic 
range by obtaining 'shallow' and 'deep' images. The field cov- 
ered by the FORS2 observations is 7'. 06 x 7'. 06 centered on 
a(2000) = ll A 08 ffl 40 i 4, 5(2000) = -60°42'51". Photometric 
calibration was achieved by using calibration plan observations 
of a Landolt (1992) standard field in the SA 100 area containing 
star #269, also reduced by the FORS2 pipeline and included in 
the data package. Extinction and colour coefficients were taken 
from the FORS2 Quality Control Web pages maintained by 
ESO (http : //www . eso . org/observing/df o/quality/). 

Instrumental photometry was carried out under IRAF by 
means of dedicated scripts that make use of tasks in the 
DAOPHOT package (Stetson 1987). Sources were first de- 
tected with DAOFIND, using first the images of isolated, un- 
saturated bright stars to determine an approximate point-spread 
function (PSF), needed for the identification of point sources. 
Photometry of all the stars detected by DAOFIND in each 
frame was then performed by defining an undersized aperture 
having as a radius the full-width at half maximum of the PSF 

1 IRAF is distributed by NOAO, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under 
contract to the National Science Foundation 
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of the reference star, and measuring the flux inside it at the 
position of each detected object using the PHOT task under 
DAOPHOT. The rest of the flux in the PSF of each star was then 
added by fitting a circularly symmetric radial profile to the dis- 
tribution of counts of each star outside the aperture, rather than 
by using a constant aperture correction. This procedure allowed 
us to remove the photometric contamination due to other stars 
located on the wings of the PSF, and to adjust to the mild varia- 
tions of the image quality across the combined images. For the 
FORS2 observations, separate catalogs of magnitudes for each 
filter were produced from the deep and shallow images, and 
the magnitudes were then compared to determine the onset of 
the non-linearity regime of the detector. The magnitude of each 
star was taken from the measurement in the deep image unless 
its magnitude was less than 0.5 mag fainter than the start of the 
non-linearity regime. For brightest stars the measurement on 
the shallow image was preferred instead. 

A merged catalog containing the results from the EMMI 
and FORS2 observations was finally prepared. The magnitudes 
measured with FORS2 were preferred when available due to 
their higher quality, and the EMMI magnitudes were used in 
the area not covered by the VLT observations. 

2.2. Spectroscopy 

The spectroscopic observations were performed in Service 
Mode with the FORS1 instrument on the VLT using multi- 
object spectroscopy (MOS) on the night of 3/4 July 2003. A 
selection of targets was made from our EMMI observations 
(which were the only ones available at the time) based on mag- 
nitudes brighter than V - 16.2 and colours bluer than Q 
-0.9. The MOS masks in FORS1 are defined by a set of mov- 
able blades, which constrains the number of slits that can be 
placed on a small area of the sky such as that subtended by 
the WRA 75 1 cluster. With this setup we could obtain spectra 
of five OB star candidates (located within V from WRA 751 
itself) in a single 600s exposure. Spectra were taken through 
the GRIS-600B+12 grism covering the spectral range between 
3450 A and 5900 A, where the diagnostic H, He, C and N lines 
used in the MK spectral type classification scheme lie. The 
GRIS-600B+12 grism is characterised by a dispersion of 1.2 
A/pix and a FWHM resolution of 4.8 A. 
The data were reduced with the standard routines in the IRAF 
APEXTRACT and ONEDSPEC packages. 

3. Identifying blue stars 

At a galactic longitude I 290°. 7 our field probes a line of sight 
providing a nearly end-on view of the Sagittarius-Carina arm, 
which extends between ~ 1 .5 kpc and ~ 8 kpc in that direction 
(e.g. Georgelin & Georgelin 1976). At the low galactic latitude 
of the field, b - -0°.4, the far end of the arm lies less than 60 pc 
from the galactic plane. It is thus expected to find in our field 
stars with highly varying degrees of extinction, leading us to 
adopt reddening-free selection criteria to identify the earliest- 
type stars in the region. 

Figure 1 shows the (B - V), (U - B) diagram of all the 
non-saturated stars in the field with colours determined to bet- 
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Fig. 1. UBV colour-colour diagram showing the positions of the non- 
saturated stars in our field having cr(B - V), cr(U- B) < 0.2. The solid 
line is the locus of unreddened main sequence stars between spectral 
types O and M5 (Bessell 1990). The dashed line is a reddening vector 
whose length corresponds to a visible extinction A v = 5 mag. 



ter than cr(B - V), cr(U - B) = 0.2. The unreddened main 
sequence locus from Bessell (1990) is plotted, as well as a 
reddening vector showing the displacement of a B3V star ob- 
scured by Av = 5 mag with a reddening law characterized by 
A v = 3.1E[B - V], E[U - B]/E[B - V] = 0.80. The adopted 
reddening ratio E[U-B]/E[B-V], also preferred by Massey & 
Thompson (1991) in their study of Cygnus OB2, is somewhat 
larger than the standard 0.72. However, it is in better agreee- 
ment with the distribution of points in the (B - V),(U - B) 
diagram, since a shallower reddening vector with a slope of 
0.72 would leave many more stars in a region inaccessible by 
a reddened star with normal intrinsic colours. We note in any 
case that the precise choice of the value of E[U - B]/E[B - V] 
has virtually no impact on our results concerning the detection 
of early-type stars. The total-to-selective extinction ratio that 
we have adopted, R v = A V /E[B - V] = 3.1, is the classical 
one and very close to the one used by Massey & Thompson 
(1991), R v = 3.0, for Cygnus OB2. The latter is in turn a 
roundoff of the value obtained by Torres-Dodgen et al. (1991), 
Rv = 3.04 + 0.09, which is virtually indistinguishable from the 
classical value given the uncertainty. 

The bulk of stars in the field have positions in the colour- 
colour diagram suggesting intrinsic colours in the range of the 
kink of the solid curve (corresponding to spectral types be- 
tween early A and early G), reddened by A v ~ 2 mag. Their 
magnitudes are consistent with their location at typical dis- 
tances of ~ 2 kpc, in the near end of the Sagittarius-Carina 
arm. A tail of redder stars, most probably cool luminous gi- 
ants at varying distances and reddenings, is also clearly visible. 
However, the most remarkable feature related to our project is 
the loose clump of stars roughly centered on (B - V) - 1.7, 
(U - B) =! 0.3, a position that is accessible only to red- 
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Fig. 2. Color-magnitude diagram showing the positions of all the non- 
saturated stars in our field. The dashed line is a reddening vector 
whose length corresponds to the displacement of a star obscured by 
Ay = 5 mag. Stars with intrinsically blue colours having Q < -0.7 are 
noted as filled circles, and the proposed WRA 75 1 cluster members 
listed in Table 1 are plotted as filled triangles. 



dened stars of the earliest types. Their location can also be 
characterized by the classical reddening-free parameter Q = 
(U - B) - 0.80(5 - V) (Johnson & Morgan 1953) 2 , where the 
stars in this clump are characterized by Q < -0.7, roughly cor- 
responding to a main-sequence spectral type B3 or earlier. 

Figure 2 shows the position of the stars in the (B - V), V 
diagram, with the stars having Q < -0.7 marked with special 
symbols. The cr(B - V), cr(U - B) < 0.2 condition prevents 
us from identifying fainter stars that may fulfill the (2-criterion, 
and the cutoff at faint V magnitudes is thus most likely not real. 

4. The WRA 751 cluster 

4.1. Imaging 

In our initial EMMI observations we identified 20 stars fulfill- 
ing the Q < -0.7 criterion. Their spatial distribution clearly 
revealed their clustering around WRA 751, as 13 of them were 
found within a radius of only 76" from WRA 751, i.e., less 
than 3% of the field imaged with EMMI. 

The deeper FORS2 observations expand this initial census 
up to the 24 blue stars in the surroundings of WRA 75 1 that we 
list in Table 1. Their positions are plotted in Figure 3. These 
stars, marked as triangles in Figure 2, are distributed around 
B - V - 1.67, which corresponds to an average extinction 
Ay ^ 6.1 towards the cluster. Their relatively narrow range of 



2 Our definition follows the usual convention in current litera- 
ture, which reverses the sign of the original definition by Johnson & 
Morgan (1953); see e.g. Massey & Thompson (1991). The coefficient 
multiplying (B - V) is E[U - B]/E[B - V] for which, as explained 
above, we prefer to adopt a value of 0.80 rather than the classical 0.72. 
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Fig. 4. The V vs. Q diagram of the intrinsically bluest stars in the 
observed fields (see text for definitions). The evolutionary tracks com- 
puted by LeJeune & Schaerer (2001, Case e with high mass loss and 
no rotation) for initial stellar masses of 9, 12, 15, 20 and 25 M Q are 
shown with dashed lines. The tracks for M; = 40 and 60 M Q are the 
dotted and dot-dashed lines, respectively. The additional long-dashed 
line represents the evolutionary track for M, = 120 M . All tracks 
have been scaled to a distance of 6 kpc and reddened by Ay = 6.1 
(see Sect. 4.1). The Zero- Age Main-Sequence (ZAMS) is traced with 
a solid line. 



colours indicates that the extinction within the cluster is rather 
small, probably not exceeding A v . ml ~ 1.5 mag. Most of the 
blue objects listed in Table 1 and plotted in Figure 3 are pro- 
jected within a distance of 1 17" of WRA 751, and the radius of 
the cluster is approximately 1 ' . Four additional stars are seen 
at larger distances, up to almost 4' from WRA 75 1 . The fairly 
strong degree of concentration towards WRA 75 1 of the other 
members suggests that these stars may not be associated to the 
cluster. However, both their colours and their magnitudes are 
within the range defined by the early stars closer to WRA 75 1 , 
and we therefore consider them as possibly related to it, also 
given the rarity of stars with such photometric characteristics 
in the rest of the field covered by our observations. We have 
thus included them in our census listed in Table 1 . 

The observed V vs. Q diagram provides an approximation 
to the H-R diagram of the cluster, and is presented in Figure 4. 
Evolutionary tracks for stars of different initial mass as given 
by LeJeune & Schaerer (2001) in Case e (with high mass loss 
and no rotation, see Meynet et al. 1994) are also plotted. Tracks 
corresponding to initial masses of 9, 12, 20 and 25 M Q are 
traced with a dashed line, while those for M, = 40 and 60 M G 
are shown with a dotted and a dot-dashed line, respectively. The 
additional long-dashed line represents the evolutionary track 
for M, = 120 M . All the tracks have been scaled to a distance 
of 6 kpc (see Sect. 5) and reddened by Ay = 6.1 (see Sect. 
4.1). Their corresponding Zero-Age Main-Sequence (ZAMS) 
is traced with a solid line. The large errors in the horizontal 
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Table 1. Likely members of the WRA 751 cluster 



5 



ID 


a 


(2000) 


5(2000) 


V 


(U 


- 


B) 




(B 


- V) 


Q 


1 


1 1 


08:37.6 


-60:42:46 


14.54 


±0.01 


0.31 


± 


0.04 


1 


59 


± 0.00 


-0.96 ± 0.04 
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± 0.05 


a m i a t a 
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J 


1 1 
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n 1 1 
0. 13 
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n nc 
0.05 
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i n n 1 
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1 aq i n nc 
— 1.U3 ± U.05 
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11 


08:41.2 


-60:42:34 


16.29 
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± 
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A Ol i A 1 T 
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11 


no . a 1 
08:41.6 


t~r\. Al . 1 o 
-60:43:18 


16.85 


±0.01 


n on 
0.80 


± 
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± 
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9 


11 


08:41.7 
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1 
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11 
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± 
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17 
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direction, mainly due to the faintness of the stars in U, pre- 
vent us from identifying accurate evolutionary phases and ini- 
tial masses, and only spectroscopy can remove this degener- 
acy. Nevertheless, Figure 4 indicates that the cluster stars can 
be quite massive, spaning a range between 9 and 120 M (the 
upper limit of the evolutionary tracks), where the lower limit of 
9 M is due to the adopted constraint of Q < -0.7. 

4.2. Spectroscopy 

The spectra of the five OB candidates observed with FORS 1 in 
July 2003 are shown in Figure 5, where the stellar continuum 
has been normalized to unity. The stars ID number is taken 
from Table 1 . 

A number of absorption lines typical of the O and B 
spectral types can be easily identified: Hy, Hel /14471, Hell 
,L14200,4686, CIII /U4070,4650 and Oil ,14367. The feature in 
absorption at /14428 is due to a diffuse interstellar band (DIB); 
its strength does not significantly vary among the observed 
stars, thus indicating quite similar reddenings. 

The spectra were classified independently by each au- 
thor via comparison with the spectral atlases of Walborn & 
Fitzpatrick (1990) and Hanson (2003). The mean spectral types 
derived for the targets are listed in Table 2, ranging from as- 
early-as 08 to as-late-as Bl. We estimate our classification 
to be accurate at less than 1.5 subtypes and two luminos- 
ity classes. We used the intrinsic UBV colour terms given by 
FitzGerald (1970) as a function of spectral type and luminosity 
class, and derived from the observed photometry of the targets 
their extinction (in Table 2). This turns out to be in quite good 
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Fig. 5. The FORS1 spectra of five OB stars in the WRA 751 cluster. 
Their stellar continuum is normalized to unity. The most relevant lines 
for spectral classification are shown together with the diffuse interstel- 
lar band in absorption (DIB) atA^ 4430 A. 



agreement with the extinction estimated from the photometry 
alone (see Sect. 4.1). 
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Table 2. The spectral classification of five early-type stars in the WRA 751 cluster, their temperature and absolute magnitude. A 
distance of 6 kpc is assumed 



ID 


Spectral 
Type 


A v 


log(T e// [K]) 
(Slesnick et al. '02) 


log(T c// [K]) 
(Massey et al. '05) 




1 


BOI 


5.67 


4.37 


4.47 


-5.02 


5 


09V 


5.49 


4.60 


4.50 


-4.40 


8 


Blla 


6.57 


4.25 




-5.61 


20 


091 


5.98 


4.49 


4.50 


-3.77 


22 


Q8V 


6.54 


4.60 


4.53 


-4.22 



The dereddened observed (B - V) colours were then used 
to determine the effective temperature (T e ff, in Table 2) of the 
targets via the empirical calibrations of Slesnick et al. (2002). 
Although these calibrations are not very sensitive to T e ff of O 
and early B stars, and do not take into account the effect of 
stellar wind blanketing, they extend to the B1I spectral type 
which is not included in the most up-to-date temperature scales 
published by Repolust et al. (2004) and Massey et al. (2005). 
Nevertheless, we have reported in Table 2 the T e ff values given 
by Massey et al. (2005) for stars earlier than BOI. The differ- 
ence between these two sets of effective temperatures is less 
than or comparable to the lcr errors on our spectral classifica- 
tion (see Sect. 5). For this reason and for the sake of homogene- 
ity, we adopt the effective temperatures given by the empirical 
calibrations of Slesnick et al. throughout Sect. 5 and in Figure 
6. 

5. Discussion 

Our observations reveal the presence of a cluster of blue stars 
spatially coincident with WRA 75 1 . Given the extreme paucity 
of early-type stars (with Q < -0.7) in the observed field, we 
identify this cluster as the birth-site of WRA 75 1 . A chance 
alignment of a cluster of massive stars with a LBV is extremely 
unlikely, given that both classes of objects are rare. 

Various distance estimates can be found in literature for 
WRA 751. van Genderen et al. (1992) derived a distance larger 
than 4-5 kpc from the photometry of field stars around 
WRA 75 1 . This is consistent with the more recent determina- 
tion of 6 + 1 kpc (Nota, private communication) based on the 
radial velocity of WRA 75 1 and its surrounding HII region. We 
have adopted the latter value in this paper; we prefer it over a 
distance simply derived from the spectral types in Table 2 and 
the available photometry in Table 1, since it is well known that 
the absolute magnitude of early-type stars displays a large scat- 
ter within any given spectral type (Garmany 1990, Jaschek & 
Gomez 1998). At this distance, the angular extent of the cluster 
translates into a radius of 3.4 pc, with WRA 751 close to the 
cluster centre. Assuming that all stars more massive than 9 M Q 
have been detected (cf. Figure 4) and that the initial mass func- 
tion of the cluster is represented by a Miller & Scalo's (1979) 
law, we roughly estimate the cluster mass to be 2.2 x 10 3 M , 
i.e. about one order of magnitude larger than the Trapezium 
cluster in Orion (Herbig & Terndrup 1986). 

The H-R diagram of the observed five OB stars in the 
WRA 751 cluster is plotted in Figure 6. The lcr errors on 



T e ff simply reflect the uncertainty of our spectral classification, 
while those on My take also into account the uncertainty on the 
adopted distance, the photometric errors and the uncertainty 
in R v (see Sect. 3). The evolutionary tracks (dashed lines) are 
taken from LeJeune & Schaerer (2001, Case e with high mass 
loss and no rotation, see also Meynet et al. 1994); their ZAMS 
is represented with a solid line. We have added to the observed 
OB stars WRA 75 1 itself, for which we have assumed a mean, 
observed magnitude V = 12.14, an extinction Ay = 5.58 mag 
and an 09.51 spectral type (cf. Hu et al. 1990). The effective 
temperature of WRA 75 1 has been derived once again from the 
empirical calibrations of Slesnick et al. 

The position of the targets in the H-R diagram indicates 
that the observed OB stars span a mass range between 12 
and 40 M Q , although the uncertainty on our spectral classi- 
fication allows masses larger than 40 M G in the case of the 
two 08V and 09V stars (# 22 and 5 respectively). The pro- 
genitor of WRA 751 may have been as massive as ~ 80 M . 
Since the brighter and bluer stars were selected for the spectro- 
scopic follow-up, we are confident that no star earlier than 08 
is present in the cluster. If the absence of stars earlier than 08 
were due to their evolution off the main sequence, this would 
set a lower limit on the age of the cluster of about 4 Myr, which 
turns out to be consistent with the current evolutionary phase of 
WRA 75 1 . We note however that the absence of such early-type 
stars may also be simply due to an insufficient sampling of the 
upper end of the main sequence (i.e., small-number statistics) 
given the limited mass that we estimate for the cluster. Thus, 
a younger age cannot be ruled out, although the classification 
of three late O/early B stars with luminosity class I argues for 
significant evolution implying an age of a few million years. 

The data presented here offer a first glimpse of the nature 
of the birth-cluster of WRA 75 1 , and need to be complemented 
with additional observations, leading to an accurate determi- 
nation of its distance, stellar content and age. Its identification 
and study are an important step towards our understanding of 
the evolution of WRA 75 1 into the Luminous Blue Variable 
phase. 
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Fig. 6. The H-R diagram of the newly discovered OB stars around 
WRA 751, at an assumed distance of 6 kpc (filled circles). The LBV 
star WRA 751 is shown with an open star symbol. Evolutionary tracks 
are from LeJeune & Schaerer (2001, Case e with high mass loss and 
no rotation). Tracks corresponding to initial masses of 12, 15, 20 and 
25 M are traced as dashed lines, while those referring to M, = 40, 
60 and 85 M Q are shown with a dotted, dot-dashed and solid line, 
respectively. 
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Fig. 3. A 6-band image of the field around WRA 751, with the stars having Q < -0.7 marked (see Sect. 3). The field covers 7'. 06 x 7'. 06 and 
WRA 75 1 is at its center. North is up and East to the left. 



